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Characterizations  of  minerogenic  particles  in  support  of  modeling  light  scattering  in 
Lake  Superior  through  a  two-component  approach 

Feng  Peng,3-*  Steve  W.  Effler,3  David  O’Donnell,3  Alan  D.  Weidemann,b  and  Martin  T.  Auerc 

11  Upstate  Freshwater  Institute,  Syracuse,  New  York 

b Naval  Research  Laboratory,  Stennis  Space  Center,  Mississippi 

c  Department  of  Civil  and  Environmental  Engineering,  Michigan  Technological  University,  Houghton,  Michigan 
Abstract 

The  role  of  suspended  minerogenic  particles  in  light  scattering  in  eastern  Lake  Superior  and  Keweenaw  Hay  (1 1 
sites)  during  July  of  2006  was  evaluated  with  an  individual  particle  analysis  technique  (scanning  electron 
microscopy  interfaced  with  automated  image  and  X-ray  analyses,  SAX),  along  with  bulk  measuren  mts  ol 
particulate  scattering  and  backscattering  coefficients  ( bp  and  bbp)  and  chlorophyll  a  concentration  ([Chi  ).  SAX 
measurements  provided  information  on  light-scattering  attributes  of  minerogenic  particles,  including  chemical 
composition  and  particle  size  distribution  (PSD).  The  data  were  used  in  Mie  theory  calculations  for  estima  ions  of 
minerogenic  scattering  and  backscattering  coefficients  ( bm  and  6bm).  The  bay  had  higher  concentral  ons  of 
minerogenic  particles  and  higher  values  of  6m,  6b<m,  6bp,  and  the  backscattering  ratio  (bhp:bp)  than  the  pelagic 
sites.  Minerogenic  scattering  was  primarily  attributable  to  clay  mineral  particles  in  the  size  range  of  1-10  /  m  The 
PSDs  deviated  from  the  (Junge)  pattern  of  monotonic  increase  of  particle  numbers  with  decreasing  parti  le  size 
The  estimates  of  bm  and  /?b>m  combined  with  those  of  organic  particulate  components,  as  represented  by  ei  lpirical 
Case  1  bio-optical  models,  showed  reasonably  good  closure  with  bp  and  6bp,  supporting  the  credibility  of  the 
SAX-Mie  approach.  Variations  in  the  bbp:bp  ratio  are  strongly  related  to  spatial  differences  in  the  dative 
contributions  of  phytoplankton  vs.  minerogenic  particles  to  scattering,  with  higher  values  observed  wl  ere  the 
minerogenic  component  is  more  important.  SAX  can  advance  the  understanding  of  optical  variab  lity  by 
establishing  the  contributions  of  well-defined  minerogenic  constituents  that  strongly  influence  optical  pre  Gerties 


Light  scattering  by  particles  is  a  fundamental  feature  of 
radiative  transfer  in  water  and  the  attendant  signal 
available  for  remote  sensing.  Two  important  metrics  of 
light  scattering  are  the  total  particulate  scattering  and 
backscattering  coefficients  (bp  and  bbp,  m-1)  that  represent 
integrations  over  the  angular  range  of  0  —  n  and  nJ2  —  n, 
respectively,  of  the  volume-scattering  function  (VSF,  /?). 
The  backscattering  ratio,  bbp:bp  (£bp),  has  recently  been 
used  as  a  diagnostic  to  infer  particle  composition  (Boss  et 
al.  2004;  Loisel  et  al.  2007,  Whitmire  et  al.  2007),  model  the 
underwater  light  field  (Mobley  et  al.  2002),  and  explain 
remote-sensing  variability  in  coastal  areas  (Lubac  and 
Loisel  2007).  These  optical  properties  depend  on  various 
attributes  of  the  particle  population,  including  the  compo¬ 
sition,  particle  size  distribution  (PSD),  particle  shape  and 
inhomogeneity,  and  number  concentration  (AO,  with  the 
exception  that  bbp  is  independent  of  N. 

Aquatic  ecosystems  represent  complex  media  with 
respect  to  particulate  scattering  because  of  the  great 
heterogeneity  of  particle  characteristics  that  influence 
related  optical  properties  differently  (Stramski  et  al. 
2004r/,  2007;  Peng  and  Effier  2007).  The  extent  to  which 
differences  in  bulk  optical  properties  can  be  explained  by 
bulk  properties  of  a  particle  population  such  as  concentra¬ 
tions  of  chlorophyll  a  ([Chi])  or  suspended  particulate 
matter  is  limited  (Babin  et  al.  2003;  Stramski  et  al.  2007). 
Moreover,  such  metrics  provide  little  insight  into  the  roles 
of  the  various  particle  types  in  regulating  these  optical 
properties. 


*  Corresponding  author:  fpeng@upstatefreshwater.org 


Stramski  et  al.  (2001,  2004r/,  2007)  h  ive  advocated  a 
reductionist  approach  of  increased  partiti  ming  of  inherent 
optical  properties  (IOPs)  into  a  ‘cei  l  i in  manageable 
number'  of  important  functional  constitn  *nts.  Though  the 
partitioning  of  bp  (or  bbp)  into  organic  and  minerogenic 
scattering  coefficients  (bQ  and  bm) 

bp  =  b0  4*  bm  ( 1 ) 

is  more  rudimentary  than  the  goals  of  the  reductionist 
approach,  such  an  effort  has  rarely  been  r  gorously  applied 
to  real  systems  (Peng  et  al.  2007;  Snyder  t  al.  2008;  Stavn 
and  Richter  2008).  A  reductionist  version  of  Eq  1  would 
further  partition  organic  and  minerogenic  particulate 
components  into  more  constituents.  Stramski  et  al.  (2001) 
implemented  the  reductionist  approach  in  a  modeling 
analysis  of  IOPs,  focusing  primarily  oi  partitioning  of 
planktonic  contributions  while  including  the  effects  of 
minerogenic  particles  through  reasonable  assumptions  of 
their  characteristics.  These  model  sim  jlations  yielded 
insights  on  optical  variability  of  the  ocean  and  suggested 
that  minerogenic  particles  were  potenth  lly  important  in 
regulating  /?bp,  occasionally  even  in  open  ocean  waters. 
Subsequently,  these  investigators  have  1  urther  evaluated 
the  potential  origins,  optical  properties,  md  implications 
for  remote  sensing  of  terrigenous  minerogenic  particles  in 
seawater  (Stramski  et  al.  2004/?,  200'  Wozniak  and 
Stramski  2004)  In  a  field  study  to  undersi  and  the  observed 
IOP  variability  in  New  England  continental  shelf  waters. 

Green  et  al.  (2003)  conducted  direct  individual  particle 

analyses  and  included  a  mineral  group  in  1  heir  five  particles 
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types.  Their  findings  emphasized  the  necessity  to  measure 
mineral  particles  for  full  understanding  of  backscattering 
variations  in  the  ocean. 

Minerogenic  particles  play  even  greater  (sometimes 
dominant)  roles  in  regulating  light  scattering  in  coastal 
Case  2  (Stramski  et  al.  2001;  Loisel  et  al  2007)  and  inland 
waters  (Kirk  1994;  Peng  and  Effler  2005;  Peng  et  al,  2007). 
In  even  the  most  dilute  and  highest  clarity  lacustrine 
systems,  minerogenic  particles  strongly  influence  apparent 
optical  properties  (AOPs;  Swift  et  al.  2006),  including 
remote-sensing  signals  (Li  et  al.  2004).  Bukata  et  al.  (1978) 
recognized  the  importance  of  suspended  minerals  in 
shaping  the  underwater  and  above-water  light  fields  of 
Lake  Superior,  and  included  an  inorganic  particulate 
component  in  the  development  of  a  semi-empirical  model. 
Recently,  an  individual  particle  analysis  (IPA)  technique, 
scanning  electron  microscopy  interfaced  with  automated 
image  and  X-ray  analyses  (SAX),  has  been  used  to  directly 
measure  the  light-scattering  attributes  (A,  PSD,  composi¬ 
tion,  and  shape)  of  minerogenic  particles  of  sizes  >0.2  pm 
in  several  inland  waters  (Peng  and  Effler  2007;  Peng  et  al. 
2007).  SAX  results  served  as  inputs  to  Mie-theory 
calculations  of  the  scattering  and  backscattering  efficiency 
factors  (gbm  and  (Jhbm,  dimensionless)  of  individual 
particles,  supporting  direct  estimates  of  spectral  bm  and 
the  minerogenic  backscattering  coefficient  ( bbftn ),  as  de¬ 
scribed  for  bm  below 


where  X  is  the  wavelength  of  light,  Nm  is  the  number  of 
minerogenic  particles  per  unit  volume  of  water  (number 
m  3),  PAm ,  is  the  projected  area  of  minerogenic  particle  /, 
and  V  is  the  sample  volume.  Likewise,  bbm(X)  was 
calculated  with  the  backscattering  efficiency  factor,  <2bbnv- 
These  two  efficiency  factors  of  a  particle  depend  on  its  size 
(d)  and  (relative)  complex  refractive  index  (w,  depends  on 
the  composition  of  the  particle),  and  on  X.  Moreover,  SAX 
supports  partitioning  of  bm  and  bbm  into  multiple 
minerogenic  particles  types,  such  as  clay  minerals,  quartz, 
and  calcite  (Peng  et  al.  2007),  consistent  with  the  goals  of 
the  reductionist  approach.  The  integrated  SAX-Mie 
method  has  been  described  as  precursor  work  for 
advancing  the  reductionist  approach  (Stramski  et  al.  2007). 

Establishment  of  the  representativeness  of  the  SAX-Mie 
estimates  of  bm  and  bbr m  is  central  to  advancing  the 
implementation  of  the  reductionist  approach.  Some  pro¬ 
gress  has  been  reported  recently  in  this  regard  for  a  limited 
number  of  particle-rich  inland  waters  (Peng  and  Effler 
2007;  Peng  et  al.  2007).  Demonstrations  for  most  natural 
systems  are  complicated  by  noteworthy  contributions  of 
organic  particles  that  are  presently  not  readily  amenable  to 
the  same  level  of  independent  individual  particle  charac¬ 
terizations.  This  issue  was  largely  avoided  in  the  study  of 
one  reservoir  (Peng  and  Effler  2007)  where  minerogenic 
particles  dominate,  and  reasonably  good  consistency 
between  SAX  Mie  estimates  of  bm  and  bulk  measurements 
of  the  attenuation  coefficient  at  660  nm,  c(660)  (a  reliable 
surrogate  of  b\  Babin  et  al.  2003)  was  reported.  In  addition, 


summations  of  SAX-Mie-based  estimates  of  bm( 650)  and 
organic  particulate  scattering,  as  estimated  by  an  empirical 
bio-optical  [Chl]-based  model  developed  for  Case  1  waters 
(Loisel  and  Morel  1998),  showed  good  closure  with  bulk 
measurements  of  />p(650)  for  several  inland  waters  enriched 
in  particles  (Peng  et  al  2007).  Similar  closure  efforts  have 
to  date  not  been  attempted  in  evaluating  the  representa¬ 
tiveness  of  SAX-Mie  estimates  of  bbr m(x).  Success  in  such 
closure  efforts  would  result  in  a  simple  two-component 
model  that  could  estimate  bp  and  bhp  and  resolve 
contributions  of  minerogenic  and  organic  particles  based 
on  SAX  and  [Chi]  measurements. 

This  paper  documents  the  light-scattering  features  of 
minerogenic  particles  assessed  by  SAX,  and  bulk  measure¬ 
ments  of  bp,  bb p,  and  [Chi],  for  pelagic  near-surface  waters 
of  oligotrophic  Lake  Superior  and  one  of  its  major 
embayments.  The  SAX-Mie-based  estimates  of  bm  and 
bbrm  are  tested  for  consistency  through  the  use  of  Case  1 
empirical  bio-optical  models  to  independent!  estimate  b0 
and  the  organic  backscattering  coefficient  (hbt  t)  from  [Chi] 
and  by  comparing  the  summations  of  the  min.irogenic  and 
organic  components  with  bulk  measurements  <  f  bp  and  bbp. 
Contributions  to  bm  and  bb  m  from  particle  ck  sses  that  are 
based  on  composition  and  size  are  resolved,  and  the 
dependence  of  spatial  differences  of  /?hp  oi  metrics  of 
composition  of  the  particulate  assemblage  s  evaluated. 
Lastly,  the  role  of  minerogenic  particles  in  influencing  light¬ 
scattering  levels  in  Lake  Superior  and  related  implications 
for  remote  sensing  are  considered.  This  paj  er  advances 
earlier  SAX-based  light-scattering  studies  (Pei  g  and  Effler 
2007;  Peng  et  al.  2007)  by  extending  to  more  dilute 
conditions,  addressing  closure  for  estimates  of  hbp  and 
demonstrating  the  interplay  between  minerot  enie  particle 
contributions  and  bbp. 

Methods 

Study  system,  sampling,  and  field  measure  nents  Lake 
Superior  is  the  westernmost  (48°00'N,  88°00  W),  largest, 
and  most  oligotrophic  of  the  Lauren tian  Gre  it  Lakes.  Its 
residence  time  is  —200  yr.  Water  clarity  is  rc  latively  high 
(Secchi  depth  15-20  m)  in  pelagic  areas  of  t  le  lake.  The 
role  of  the  minerogenic  particle  populat  on  in  light 
scattering  in  the  lake  has  not  been  characterized.  Lake 
Superior  has  been  described  as  optically  com  >lex  (Case  2; 
Li  et  al.  2004).  In  particular,  absorption  by  c  iromophorie 
dissolved  organic  matter  (CDOM)  is  high  rel  itive  to  that 
associated  with  phytoplankton  pigments  (Gors  et  al  2008) 
and  sediment  is  believed  to  contribute  to  1  >cal  remote¬ 
sensing  signatures  (Budd  2004),  though  lake  wide  effects 
are  uncertain  (Li  et  al.  2004). 

A  total  of  1 1  sites  (Fig.  1)  were  visited  ove  the  interval 
of  05-14  July  in  2006.  These  included  two  lot  itions  in  the 
Keweenaw  Waterway  (KW)  proximate  to  a  iru  jor  tributary 
inflow  (the  Sturgeon  River;  Fig.  1),  six  sites  along  the  axis 
of  Keweenaw  Bay  (KB),  and  three  pelagic  sii  *s  in  eastern 
Lake  Superior  (ELS).  Site  KW1  reflects  the  k  :al  effects  oi 
the  Sturgeon  River  inflow  to  Keweenaw  Bay  through  the 

Keweenaw  Waterway.  Site  KW2  is  a  proximate  (within 
1  km)  nearshore  bay  location,  included  to  pro  ide  contrast 
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Keweenaw  Waterway  (KW;  two  sites),  Keweenaw  Bay  (KB;  six 
sites),  and  Eastern  Lake  Superior  (ELS;  three  sites). 

to  KW1.  The  other  bay  sites  extended  along  a  transect  from 
the  inner-(KBl)  to  outer-bay  (KB6)  bounds  (^33  km).  The 
three  pelagic  sites  in  eastern  Lake  Superior  (ELS;  shelf  and 
profundal  locations)  extend  over  —155  km  from  Michipi- 
coten  Island  (Fig.  1).  At  each  site  in  situ  optical  measure¬ 
ments  were  conducted  and  near-surface  samples  were 
collected  for  IPA  and  [Chi]  (Parsons  et  al.  1984)  analyses. 

Spectral  attenuation,  absorption,  and  backscattering 
coefficients  [c(A),  n(A),  and  bb(X):  m_I]  were  measured  with 
a  combined  instrument  package  (single  profiling  cage, 
integrated  outputs)  consisting  of  an  ac-s  (spectral  absorp¬ 
tion  and  attenuation  meter)  and  a  BB9  (backscattering 
meter)  from  WET  Labs.  Backscattering  measurements 
were  unavailable  for  the  two  KW  sites.  The  ac-s  measures 
c(X)  and  a(X)  (25-cm  path  length)  relative  to  pure  water, 
with  a  spectral  resolution  of  4  nm  over  the  range  of  400- 
730  nm  This  instrument  was  calibrated  according  to 
guidelines  of  the  manufacturer  with  pure  water  and 
checked  between  deployments  following  the  manufacturer’s 
calibration  protocol  (WET  Labs,  Inc.  2006).  Adjustments 
for  the  effects  of  temperature  on  pure  water  absorption  and 
attenuation  were  made  according  to  Pegau  et  al.  (1997). 
The  absorption  spectra  were  corrected  for  scattering  error 
following  the  manufacturer's  protocols.  The  spectral 
particulate  scattering  coefficient,  />P(A),  was  determined  by 
difference,  c(a)  —  a(A).  The  BB9  measured  VSF  at  117° 

[/?p(l  1 7  ,  /)]  at  nine  discrete  wavelengths,  for  which  c{X) 
and  a{X)  measurements  were  made.  The  angle  of  117°  has 


been  identified  as  a  minimum  convergence  point  for 
variations  in  P(0,  X)  caused  by  suspended  particles. 
Backscattering  [/?bp(2)]  was  calculated  from  /Ip(117°,  X) 
according  to  />bp(A)  “  27Tx/?p(l  17°, A),  where  x  “  1  1  (Boss 
and  Pegau  2001).  Corrections  for  attenuation  along  the 
path  from  the  light  source  to  the  sample  volume  and  then 
to  the  detector  were  made  based  on  a(/.)  measurements 
from  the  ac-s  (WET  Labs,  Inc.  2007).  Values  of  bp  and  bb p 
reported  here  are  averages  of  those  collected  over  a  0.5-m 
depth  interval  in  near-surface  waters. 

SAX  protocols — Sample  preparation  and  instrumenta¬ 
tion  methodology  for  SAX  have  been  described  previously 
(Peng  and  Effler  2005,  2007),  and  are  onl}  briefly  reviewed 
here.  Suspended  particles  from  known  volumes  of  water 
samples  were  deposited  onto  polycarbonate  membranes 
(OA-pm  pore  size)  by  pressure  filtration,  air-dried,  and 
coated  with  carbon.  The  SAX  instrum  ntation  was  an 
Aspex  PSEM  2000  system  controlled  by  automated  feature 
analysis  (AFA)  software  (Aspex®).  Abou  2000  individual 
particles  were  analyzed  for  each  of  the  10  s  imples  (no  water 
sample  from  KW2  was  available)  for  both  morphology  and 
composition  characteristics  (~3  h  per  sail  pie). 

Based  on  a  rotating  chord  algorithm,  AFA  represents 
each  particle  by  a  series  of  radiating  triang  es  formed  by  the 
particle's  centroid  and  16  ehords  drawn  through  the 
centroid  to  the  edges  of  the  particle  at  11°  increments. 
The  projected  area  (PA,-;  m2)  of  a  particle  i:.  measured  as  the 
sum  of  all  triangular  areas.  Particle  size  (<«,;  /mi)  is  defined 
as  the  circular  area  equivalent  diameter.  Particle  shape  is 
described  here  by  the  aspect  ratio  (ASP),  the  ratio  of  the 
length  of  the  longest  chord  to  that  of  its  o  thogonal  chord. 
The  ASP  value  is  a  measure  of  a  particle* s  deviation  from 
sphericity  (ASP  of  a  sphere  is  1). 

Composition  of  individual  minerogt  nic  particles  is 
assessed  by  SAX  from  acquisition  of  X-  ay  counts  of  16 
elements  (Na  and  higher  atomic  numbers  including  Al,  Si, 
Ca,  Fe,  Mn),  and  is  represented  by  its  non  lalized  elemental 
X-ray  counts.  Based  on  their  elemental  X-ray  composition, 
minerogenic  particles  w'ere  sorted  into  five  predefined 
generic  minerogenic  particle  types  (Peng  rt  al.  2007),  with 
assigned  values  of  tn  {in  =  n  -  n' i ,  where  n  and  n'  are  the 
real  and  imaginary  parts  of  m):  ‘Clay’  (aluminum  silicates, 
m  =  1.173  -  0.00 1  /);  ‘Quartz’  (mineral  Sr  )2,  m  1.155 
0.0001/);  ‘Si-rich’  (Si-containing  minerals  assumed  m 
1.173  —  0.001/);  ‘Ca-rich’  (particles  contain  ^40%  Ca  X- 
rays,  m  =  1.19  -  0.0005/);  and  ‘Mist  (miscellaneous 
inorganic  particles,  assumed  m  =  1.173  -  0.001/). 
Specification  of  m  values  for  these  particle  types  was  based 
on  literature  compilations  (Kerr  and  Rogers  1977;  Woz- 
niak  and  Stramski  2004).  ‘Diatom’  particles  were  differen¬ 
tiated  from  ‘Quartz’  (both  have  >90%  Si  "(-rays)  based  on 
their  lower  X-ray  densities  (Peng  et  al.  2002).  Hereafter 
particles  in  the  ‘Clay’  class  are  describee  as  clay  mineral 
particles  (e.g.,  kaolinite),  a  chemical  groi  ping  that  differs 
from  the  use  of  ‘clay’  to  describe  a  partic  e  size  group. 

A  useful  bulk  property  of  an  analyzed  s  mple  is  the  total 
particle  projected  area  concentration  or  the  particle 

projected  area  per  unit  volume  of  water  (PAV;  m-1),  a 
metric  that  is  approximately  proportional  o  the  magnitude 
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of  light  attenuation  (Peng  and  Effler  2007)  and  remote¬ 
sensing  reflectance  (Sydor  and  Arnone  1997).  SAX 
measurements  support  determination  of  the  PAV  associat¬ 
ed  with  each  of  the  specified  minerogenic  particle  types. 
The  minerogenic  component  of  PAV  (PAVm)  is  the 
summation  of  values  for  the  five  specified  particle  types. 

PSDs  and  scattering  calculations —  PSDs  of  minerogenic 
particles  are  presented  as  a  size  distribution  function  (or 
density  function),  F(d)  =  1 IV  X  ANIAd ,  where  Vis  the  total 
number  of  minerogenic  particles  in  V,  and  AN  is  the 
number  of  mineral  particles  in  a  size  interval  of  Ad.  Size 
distributions  were  extended  to  particles  >0.2  pm,  a  size 
somewhat  smaller  than  the  nominal  pore  size  of  the  used 
filters.  Atteia  et  al.  (1998)  demonstrated  that  particles  much 
smaller  than  the  respective  pore  sizes  of  the  membrane 
filters  were  retained.  In  paired  measurements  of  samples 
collected  on  0.2-^m-  and  0.4-//m-pore-sized  filters,  we  have 
found  that  particles  in  the  0.2-0.4-/tm  size  range  were  only 
slightly  underrepresented  by  the  0.4-/m-pore-sized  filter, 
and  that  the  general  PSD  patterns  were  not  affected  by  the 
filter  pore  size.  Results  presented  here  for  the  0.2-0.4-/m 
size  interval  are  not  presented  as  quantitative  but  serve  to 
depict  the  general  pattern  for  these  small  particles. 

The  minerogenic  particle-scattering  and  backscattering 
coefficients  (6m  and  /?bm)  were  calculated  from  SAX  results 
according  to  Eq.  2.  Values  of  the  efficiency  factors  Qbmj 
and  (2bbm./  were  obtained  from  the  BHMie  algorithm  by 
Bohren  and  Huffman  (1983),  based  on  Mie  theory  for 
homogenous  spheres.  Scattering  calculations  and  field 
observations  are  reported  here  for  X  =  650  nm;  reference 
to  wavelength  is  omitted  hereafter. 

Adjustments  of  PA/  values  were  made  to  accommodate 
the  effect  of  the  preferred  orientation  (lying  flat)  of  mineral 
particles  on  the  filter  membranes  (Jonasz  1987),  based  on 
the  ASP  values  of  the  particle;  if  ASP  >1.5,  then 
PA*  =0.8  x  PA/,  where  PAJ  is  the  adjusted  PA  value  for 
the  /th  particle  (Peng  et  al.  2007).  Average  scattering 
efficiency  factors  of  the  minerogenic  particle  populations 
(^b.n^  and  <2bb.m>)  of  each  site  were  calculated  as  the 
ratios  of  bm  and  /?bm  to  PAVm.  These  average  efficiency 
values  reflect  the  combined  effects  of  PSD  and  chemical 
composition. 

Empirical  bio-optical  models  for  scattering  (Loisel 
and  Morel  1998)  and  backscattering  (Morel  and 
Maritorena  2001)  coefficients  as  functions  of  [Chi] 
for  Case  1  (open  ocean)  waters  were  applied  here  as 
estimators  of  the  organic  components  of  bp  and  bhp  (i.e.,  b0 
and  bh. «). 

/>o(650)  =  0.347[Chl]O766(660/650)0.97  (3) 

/>b.o(650)  =  {0.002  +  0.0 1  [0.50  —  0.251og10[Chl]]  (650/ 550)v} 
x  /}o(550)  (4) 

where 

j 0.5(log10[Chl]  —  0.3)  for  0.02 <  [Chi]  <2 mg m-3 
V_  (0  otherwise  '  a' 


and 

bo(550)  =  0.4 16  [Chi] 0  766  (4b) 

The  0.97  adjustment  in  Eq.  3  accommodates  the  effect  of 
replacing  c  in  the  original  model  with  b.  The  summations  of 
minerogenic  and  organic  scattering  and  backscattering 
components  were  compared  to  bulk  measurements  of  bp 
(Eq.  1)  and  bbp  to  evaluate  the  consistency  of  this  forward 
modeling  approach  and  the  extent  of  closure. 

Results 

Light-scattering  attributes  of  minerogenk  particles 
Light-scattering  features  of  minerogenic  particles  collected 
from  Lake  Superior,  as  characterized  by  the  SAX-Mie 
approach,  are  presented  in  Table  1.  The  co  nposition  of 
minerogenic  particles  is  represented  in  the  cot  text  of  light¬ 
scattering  implications  by  the  relative  contribution  of  the 
five  particle  types  to  PAVm.  The  ‘Misc’  typt  contributed 
<12%  to  PAVm,  indicating  that  the  adopted  classification 
scheme  performed  well  in  categorizing  the  minerogenic 
particles  in  the  lake.  Generally  uniform  com  oosition  was 
observed  (Table  1).  Clay  minerals  (‘Clay’  clas ;)  dominated 
PAVm  for  all  the  sites,  representing  an  average  of  73% 
(coefficient  of  variation,  CV  =  12%)  of  the  tc  Lai .  The  next 
largest  contributors  to  PAVm  were  the  *Quc  rtz’  and  ‘Si- 
rich’  classes  with  the  exception  of  site  ELS?  where  ‘Ca- 
rich’  was  the  second  largest  component. 

Nonspherical  morphology  prevailed  for  the  minerogenic 
particle  populations  of  all  the  sites,  as  depictec  by  the  mean 
ASP  values  (Table  1).  Though  the  shapes  of  the  minero- 
genic  particles  were  diverse  (Fig.  2A),  as  ind  cated  by  the 
relatively  large  standard  deviation  (SD)  values,  the  mean 
ASP  values  were  within  a  narrow  range  of  1.93-2.30, 
suggesting  a  similarity  in  the  distributer  of  shapes. 
However,  these  deviations  from  sphericity  ol  the  minero¬ 
genic  particles  were  modest  relative  to  those  abserved  for 
diatoms,  whose  large  mean  ASP  values  (‘ everal  times 
greater  than  those  for  the  minerogenic  particle  populations; 
Table  1)  were  attributed  mostly  to  pennate  forms  (e.g.. 
Fig.  2B). 

Representative  minerogenic  PSDs  are  prt  rented  for  a 
bay  (KB6)  and  an  open  lake  site  (ELS3;  Fig.  3  Monotonic 
increases  in  number  density  with  the  decrea  e  in  particle 
sizes,  a  pattern  commonly  described  by  the  Ju  lge  function, 
were  not  observed.  Instead,  the  PSDs  exl  ibited  right- 
skewed  patterns.  Peak  densities  were  betweei  0.4  pm  and 
0.6  pm,  around  which  greater  and  longer  tai  *d  decline  in 
particle  numbers  for  particles  larger  than  the  modal  sizes 
prevailed.  If  the  Junge  function  was  fitted  for  particles 
larger  than  the  peak  density  size,  the  slope  vah  es  (~~3)  were 
consistent  with  those  reported  in  the  literati  re  (Stramski 
and  Kiefer  1991).  Differences  among  sites  wer ;  observed  in 
number  concentrations  rather  than  the  gene  a  1  shapes  of 
PSDs  (Fig.  3).  Moreover,  similar  patterns  prevailed  for  the 
various  particle  groups,  as  illustrated  for  the  day  minerals 
at  site  KB6. 

Calculated  size  dependencies  of  bm  (Fig.  4 A)  and  bbm 

(Fig.  4B)  are  presented  for  the  same  two  sa  uples  in  the 
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Table  1  Composition,  morphology,  size  distribution,  and  scattering  characteristics  of  suspended  particle  p<.  pulations  in  Lake 
Superior  (no  observations  for  site  KB2). 


Site 

PAVm* 

(m->) 

PAVm  type  distribution  {%) 

mean  ASP±SDf 

<Qb.m> 

<2bb.ni> 

'(so.bin 

(/'Hi) 

'(50.bi.111 

(/mi) 

Clay 

Quartz 

Si-rich 

Ca-rich 

Mi  sc 

Mineral 

Diatom 

KW1 

0.182 

75.3 

7.1 

4.7 

11 

11.9 

1 .93±  1.32 

7.63±  14.5 

2.23 

0.044 

3.63 

4.17 

KW2 

0.031 

76.9 

3.1 

15.3 

1.0 

3.7 

2.04±2  10 

8.76±9.87 

2.16 

0.040 

2.08 

3.1 1 

KB1 

0.067 

74.5 

10.5 

7.0 

0.6 

7.3 

1 .991:1.53 

4.97±7.57 

2.31 

0.045 

2.23 

3.18 

KB3 

0.080 

75.0 

3.2 

18.3 

0.0 

3.6 

2.30±2.01 

1 1 .56±  13.8 

2.24 

0.051 

2.95 

3.77 

KB4 

0.030 

73.3 

11.0 

7.5 

1.2 

7.0 

1.94  ±1.44 

10.28±  17.7 

2.31 

0.047 

2.31 

3.28 

KBS 

0.043 

82.1 

3.0 

10.8 

0.1 

4.1 

2.27±2.62 

10.25±  1 1  8 

2.11 

0.038 

2.24 

3.65 

KB6 

0.050 

69.6 

15.4 

11  1 

0.0 

3.8 

1.93±  1.42 

13.12±  18.6 

2.38 

0.045 

2.05 

2.93 

ELS1 

0.016 

64.3 

5.8 

13.0 

5.7 

11.2 

1.94±2.42 

9.45±  15.0 

2.36 

0.050 

2.39 

3.39 

ELS2 

0.026 

73.3 

9.1 

8,9 

3.4 

5.3 

2  00±  1.51 

7.72±  12.8 

2.21 

0.043 

2.26 

3.33 

ELS3 

0.022 

54.7 

6.3 

14.8 

17.7 

6.5 

2.12±2.29 

10.58±  12.0 

2.38 

0.047 

2.19 

3.13 

PAVm:  total  projected  area  of  minerogenic  particles  per  unit  volume  of  water, 
f  ASP:  aspect  ratio:  SD:  standard  deviation. 


format  of  cumulative  contributions  of  particle  size  classes 
to  scattering.  Submicron  particles  and  those  larger  than 
10  pm  made  only  minor  (<10%)  contributions  to  both  bm 
and  £b<m.  The  sizes  that  correspond  to  the  50th  percentile  of 


Fig.  2.  Micrographs  of  selected  particles  by  scanning  elec¬ 
tron  microscopy:  (A)  inorganic  particles,  including  a  plate-like 
clay  mineral,  and  (B)  pennate  diatoms  (with  high  ASP  values). 
Particles  in  (A)  and  (B)  were  from  the  sample  collected  from  KW2 
on  05  July  2006. 


bm  (or  median  size  of  scattering,  d50, bm;  1  ig.  4A)  and  /?bt)11 
(J5o,bbmj  Fig-  4B)  are  useful  metrics  in  the  evaluation  of  the 
relative  importance  of  various  particle  si  ,es  to  scattering. 
These  values  were  the  highest  in  the  ini  o\v  zone  of  the 
tributary,  site  KW1  (Table  1).  Little  variat  on  was  observed 
for  the  other  sites,  with  the  values  of  rf50, bi  m  shifted  ~  1  pm 
larger  than  those  of  dso,bm  (Table  1);  th  ;  average  values 
were  3.39  pm  and  2.43  pm ,  respectively  No  systematic 
spatial  structure  in  these  metrics  was  obs<  rved 

Spatial  patterns  of  PA  Vm,  bm,  bfJ  /n,  a)  d  bulk  measures 
of  scattering — Measurements  of  PAVn  (Fig.  5A)  and 
coupled  estimates  of  bm  (Fig.  5B)  and  6b  n  (Fig.  5C)  were 
the  highest  in  the  tributary  inflow  zone  (site  KW1)  and  the 
lowest  in  the  open  waters  of  eastern  Lake  Superior  (ELS 
sites).  The  mean  PAVni  for  the  five  Keweenaw  Bay  (KB) 
sites  was  >two-fold  greater  than  for  the  three  ELS  sites. 
There  were  substantial  spatial  difference  in  minerogenic 
particle  levels  and  scattering  within  be  h  the  bay  and 
eastern  Lake  Superior.  The  influence  of  the  monitored 
tributary  input  (KW1;  Fig.  1)  on  spatial  structure  in  the 
bay  at  the  time  of  sampling  was  small  (Fig  5A).  The  spatial 


Particle  size,  d  (pn  ) 

Fig.  3.  Examples  of  mineral  PSDs  expressed  in  density 

function.  Samples  were  from  KB6  and  ELS3  sites,  and  PSD  of 
clay  minerals  (‘Clay’  class)  from  KB6  is  also  s  town. 
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Fig.  4.  Size  dependency  patterns  of  scattering  for  two  sites, 
KB6  and  ELS3:  (A)  cumulative  bm ,  and  (B)  cumulative  bbm.  The 
sizes  that  correspond  to  the  50th  percentile  of  bm  and  bbm  {clso.bm 
and  ^5o»bbm»  respectively)  are  illustrated. 

patterns  in  the  calculated  values  of  bm  and  bbm  (Fig.  5B,C) 
tracked  linearly  those  observed  for  PAVm  (Fig.  5A). 
Accordingly,  the  values  of  <(2b,m>  and  <Qbb,m>  were 
uniform;  the  average  values  were  2.23  and  0.045,  respec¬ 
tively,  with  CV  values  of  4.2%  and  9.0%.  Levels  of  [Chi] 
were  somewhat  higher  at  the  ELS  sites  than  in  the  bay  (Fig. 
5A);  average  concentrations  were  0.69  mg  m~3  and  0.63 
mg  m~3,  respectively. 

Bulk  measurements  of  bp  demonstrated  more  modest 
differences  between  the  bay  and  ELS;  the  average  for  the 
bay  (0.32  m  !)  was  only  16%  higher  (Fig.  5B).  However, 
bulk  measurements  of  bbp  (Fig.  5C)  were  on  average 
(0.004  m_1)  substantially  higher  in  the  bay  than  at  the 
ELS  sites,  to  an  extent  (two-fold)  that  approached  the 
differences  calculated  for  bbm.  Differences  in  PAVm 
explained  76%  of  the  observed  variability  in  bbp  according 
to  linear  least-squares  regression  (p  <  0.01);  however,  these 
were  not  a  significant  predictor  of  bp  dynamics. 

Performance  of  the  two-component  model-  Summations 
of  minerogenic  and  organic  scattering  contributions  are 
compared  to  the  bulk  measurements  (Fig.  5B,C).  These 
summations  were  close  to  the  measurements  of  bp\  the 
average  deviation  was  8.2%.  All  but  one  (ELS2)  of  the 
summations  of  the  independent  model  estimates  were 
within  13%  of  the  observations  (excluding  ELS2,  average 
percentage  error  of  —6%).  Though  the  linear  relationship 
between  measured  values  of  bp  and  the  summations  was 
heavily  weighted  by  the  KWl  site,  an  approach  towards 


Fig  5.  Results  of  laboratory  analyses,  in  situ  optical 
measurements,  and  light-scattering  modeling  for  '.tudy  sites  in 
Lake  Superior:  (A)  PAVm  and  [Chi],  (B)  measurer  and  modeled 
bpy  (C)  measured  and  modeled  hbp ,  (D)  backscatter  ig  ratio,  from 
measurements  and  (bbo  +  bh  m) :  (b0  +  />m),  and  (E)  lackscatlcring 
ratios  of  organic  and  minerogenic  components. 

equivalency  (slope  0.8)  was  still  indicated  for  t  le  other  sites 
(Fig.  6A,  and  inset). 

The  higher  PAVm  and  organic  scatterii  g  (Fig.  5A) 
contributed  nearly  equally  to  the  elevated  bp  measured  at 
KWl  (Fig.  5B).  Levels  of  bp  in  KB  sites  and  ELS  sites  were 
not  significantly  different  from  each  other.  Hi  wever,  there 
were  noteworthy  spatial  differences  in  the  re  ative  contri¬ 
butions  of  bQ  and  bm  to  bp  (Fig.  5B).  For  the  ;hree  pelagic 
ELS  sites,  b0  accounted  for  an  average  of  9  •%  of  bp.  In 
contrast,  its  contribution  decreased  to  a  mean  of  63%  (CV 
=  16%)  for  the  other  sites.  Conversely,  bm  \  as  of  minor 
importance  to  bp  at  the  ELS  sites  (—17%  on  iverage),  but 
contributed  more  substantially  (—34%  on  average)  at  the 
KW  and  KB  sites,  particularly  at  KB1  and  KB^  (—45%). 

These  results  support  the  use  of  the  empirica1  relationship 
(Eq.  3)  for  bv  in  this  waterbody. 
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Fig.  6.  Comparisons  of  modeling  results  (open  symbols)  with 
scattering  measurements,  with  linear  regression  analyses  (solid 
line:  linear  fit,  dashed  line  95%  confidence  limit):  (A)  bp 
(excluding  KWI  point),  with  inset  for  all  data  points  (y  =  1.10.x 
0.04,  R2  =  0.98,  p  <  0.001),  (B)  6bp,  and  (C)  backscattering 
ratio.  Filled  symbols  in  (B)  represents  summation  result  with  bbo 
predicted  by  an  alternative  model  (Huot  et  al  2008),  and  the 
corresponding  linear  regression  (lines  not  shown)  is  y  =  0.87x  + 
0  0005  ( R 2  =  0.72,  p  <  0.01). 

On  average,  the  two-component  model  estimates  of  bhp 
(b b.o  +  bhm)  exceeded  bulk  measurements  by  14% 
(Fig.  5C).  The  overeslimation  was  more  pronounced  for 
ELS  sites  than  for  the  five  KB  sites,  by  an  average  of  63% 
vs.  10%.  This  overestimation  was  mainly  caused  by  two 


outliers  in  ELS1  and  ELS2,  for  which  the  [Chl]-based 
estimates  of  bbo  alone  were  greater  than  measured  bbp  by 
18%  and  30%  (Fig.  5C),  respectively.  Nonetheless,  the  two- 
component  model  performed  reasonably  well  in  represent¬ 
ing  the  bulk  measurements  of  bbp  (Fig.  6B).  The  predictions 
explained  68%  (p  =  0.01)  of  the  variability  in  the 
observations,  according  to  linear  least-squares  regression, 
with  general  tendency  toward  equivalency  (slope  value  of 
0.9)  and  95%  confidence  limits  nearly  encompassing  all  of 
the  data. 

The  relative  importance  of  the  minerogenic  component 
was  shifted  higher  for  6bp  than  for  bp\  on  average,  hbm 
accounted  for  42%  of  the  bulk  measurements  in  the  ELS 
sites  vs.  64%  overall  for  KB  sites  (exceeded  75%  of  the  hbp 
values  for  the  KB1  and  KB3  locations).  The  highest 
predicted  bbp  value  (0.012  m_1)  was  for  KWI,  with  a 
minerogenic  contribution  of  67.4%  (Fig.  ;  C). 

Backscattering  ratio  and  its  dependencie  The  backscat¬ 
tering  ratio,  bb p,  based  on  in  situ  measurements,  was 
substantially  lower  at  the  ELS  sites  (mean  =  0.0075)  than 
the  KB  sites  (mean  =  0.0125;  Fig.  5D  The  estimated 
ratios,  based  on  SAX  results  and  bio-optic  d  models  [(bb  o  + 
^b,m)  •  (^o  +  &m)]>  on  average  were  12%  (  ange  =  7-26%) 
higher  than  those  based  on  bulk  measuren  ents  at  KB  sites. 
However,  deviations  were  more  pronoun  ced  for  the  ELS 
sites  (45%  on  average;  Fig.  5D),  associated  with  relatively 
large  predicted  contributions  by  bbo.  Tin  estimated  ratio 
for  KWI  was  higher  than  that  for  KWI.  The  respective 
backscattering^  ratios  of  the  organic  a  id  minerogenic 
components  (bbo  and  Ab  m)  are  presente  1  in  Fig.  5E  for 
reference  Values  of  bb  m  were  uniform  (m  ean  =  0.022,  CV 
=  7.7%)  and  on  average  more  than  two  t  mes  higher  than 
those  of  bbo  (mean  =  0.009,  CV  =  4.7%;  Fig.  5E), 
established  by  the  adopted  bio-optic; 1 1  models.  The 
modeled  bbp  approached  equivalency  (slo.)e  value  of  1.02, 
R2  =  0.75,  p  <  0.01,  intercept  not  signi  icantly  different 
from  0)  with  values  based  on  bulk  measun  ments  (Fig.  6C). 

The  dependence  of  bbp  values  on  in*  icators  of  bulk 
particle  composition,  [Chi]  and  the  [Chi]  p  ratio  (cp  is  the 
particulate  attenuation  coefficient  at  A  =  650  nm;  Boss  et 
al.  2004;  Loisel  et  al.  2007),  was  evaluated.  The  ratio  bbp 
was  negatively  correlated  ( R  =  0.81)  witl  log([Chl])  (Fig. 
7A).  Strong  negative  relationships  also  existed  between 
bb p  and  [Chi] :  cp  based  on  both  measurem  ents  (Fig.  7B;  R2 
=  0.88)  and  modeling  (Fig.  7C;  R2  =  0.66).  The 
dependence  of  bbp  (from  bulk  measurements)  on  the  ratio 
of  bm :  cp,  was  observed  to  be  a  strong  positive  relationship 
(R2  =  0.67;  Fig.  7D). 

Discussion 

Minerogenic  particles  in  Lake  Superior  Levels  of  bp  in 
the  pelagic  waters  of  Lake  Superior  and  ii  Keweenaw  Bay 
were  low  relative  to  most  of  the  lacustrine  systems  included 
in  the  listing  of  Kirk  (1994).  The  observed  levels  of  PAVm 
were  comparable  to  those  reported  for  5  kaneateles  Lake 
(one  of  the  Finger  Lakes  of  New  York),  ih  ;  system  with  die 

highest  clarity  among  those  characterize  i  previously  by 
SAX  (Peng  et  al.  2007).  Minerogenic  part  elcs,  mostly  clay 
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Fig.  7.  Evaluations  of  backscattering  ratio  dependence  (at 
650  nm):  (A)  observed  bbp  as  a  function  of  [Chi]  and  power-law 
fit  with  y}  reduced  to  2.7  X  10-6,  (B  and  C)  bbp  (observed  and 
modeled)  as  a  function  of  the  [Chi] :  cp  ratio,  and  (D)  observed  bbp 
as  a  function  of  the  hm :  cp  ratio.  Linear  least-squares  regression 
results  are  shown  in  (B),  (C),  and  (D). 


minerals  in  the  size  range  1-10  /im,  were  an  important 
component  of  backscattering  in  Lake  Superior  and 
Keweenaw  Bay  during  the  study,  and  were  a  significant 
contributor  to  bp  in  the  bay.  The  higher  levels  of  PAVni  in 
KB  compared  to  the  open  waters  of  the  ELS  sites  (Fig.  5A 
C)  presumably  reflected  local  inputs  of  mineral  particles 
from  the  tributary  (Sturgeon  River)  and  perhaps  resuspen¬ 
sion  Both  levels  and  spatial  heterogeneity  of  this  scattering 
component  were  probably  conservatively  represented  by 
the  summertime  sampling  of  this  study.  For  example, 
>70%  of  the  annual  sediment  loading  to  the  Keweenaw 
area  is  received  in  spring  (Auer  and  Gatzke  2004) 
Substantial  sediment  plumes  have  been  described  for  the 
lake,  in  spring  from  tributary  inputs  (Budd  2"04;  Stortz  et 
al.  1976).  The  documented  differences  in  PA\  m  within  the 
bay  and  ELS  sites  supports  the  position  thai  such  spatial 
structure  in  minerogenic  particle  concentrations  is  com¬ 
mon.  The  elevated  levels  of  bbp  and  bbm  in  Keweenaw  Bay 
represent  an  appropriate  target  for  resolu  ion  through 
remote-sensing  techniques. 

The  high  levels  of  bb  m  and  cp  (Loisel  and  Morel  1998) 
relative  to  the  prevailing  levels  of  [Chi]  and  their  generally 
noncovarying  character,  as  well  as  the  high  (  DOM  levels 
(Gons  et  al.  2008),  make  the  lake  and  its  ei  ibayments  a 
Case  2  system  Thus,  it  is  not  surprising  that  application  of 
bio-optical  remote-sensing  algorithms,  developed  for  Case 
1  waters,  resulted  in  overestimations  of  [Cl  1]  (Li  et  al. 
2004).  Wozniak  and  Stramski  (2004)  pn  dieted  such 
overestimates  from  the  effects  of  minerogenic  particles 
through  Mie  scattering  calculations  for  differs  nt  PSDs  and 
refractive  indices  (m).  Moreover,  these  invt  stigators  re¬ 
ported  that  PSD  and  m  must  be  taken  into  ac  mint,  rather 
than  just  the  common  bulk  gravimetric  me  isure  of  the 
minerogenic  particle  population,  to  quantify  hese  effects 
SAX  results  may  be  valuable  in  developing  a  more  robust 
mechanistic  algorithm  for  Lake  Superior  and  other  Case  2 
systems. 

Early  comparisons  to  other  systems — A  though  the 
number  of  freshwater  systems  studied  by  S  AX  remains 
modest,  the  general  shape  of  the  PSDs  (Fig.  )  has  been  a 
recurring  feature  (Peng  and  Effler  2007;  Peng  et  al.  2007). 
The  observed  PSD  pattern  is  consistent  with  results  from 
model  simulations  of  particle  aggregation  an  1  sedimenta¬ 
tion  in  a  closed  system  and  with  the  cor  elusion  that 
particles  (individuals  or  aggregates)  in  the  she  range  0.1 
10  /zm  are  the  most  stable  particulate  materi  il  in  aquatic 
systems  (Buffle  and  Leppard  1995;  Atteia  1998).  The 
dominant  role  of  particles  in  the  size  range  1  10  /zm  in 
regulating  bm  and  bb m  in  Lake  Superior  an  l  Keweenaw 
Bay  (Fig.  4)  is  also  documented  in  earlier  stud  es.  Observed 
differences  in  the  relative  importance  of  vai  ous  particle 
sizes  within  this  range,  as  represented  by  clso  \  alues,  reflect 
differences  in  the  details  of  the  various  PSDs  in  these 
studies.  The  average  d50  for  bm  for  this  large  1  ike  (2.4  /zm) 
was  about  one-half  of  those  reported  for  r  mch  smaller 
lakes  and  a  river  by  Peng  et  al.  (2007).  A  750  value  of 
~  1  //m  was  reported  for  typical  marine  mineral  particles 

(Babin  et  al,  2003)  based  on  Mie  theory  calculations 

(specifications  of  n  =  1.18,  Junge  function  si  ipe  range  = 
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3. 8-4. 2).  The  differences  among  these  systems  arc  superfi¬ 
cially  consistent  with  differences  in  sorting  (e.g.,  settling 
propensity)  expected  for  these  environments  (e.g.,  from 
large  rivers  and  small  lakes  to  large  lakes,  to  the  oceans). 
Both  PSD  and  particle  composition,  features  that  may  vary 
greatly  among  systems,  affect  the  size  dependency  patterns 
of  bp  and  bb p.  The  differences  in  the  cumulative  percent 
contribution  per  size  class  to  bm  and  7>bjm  reported  here 
(Fig.  4;  see  also  d$o  values  in  Table  1)  are  manifestations  of 
the  high  refractive  indices  and  PSDs  of  the  minerogenic 
particle  population  of  Lake  Superior. 

Clay  minerals  have  been  a  prominent  component  of 
minerogenic  scattering  in  all  the  characterized  inland 
waters.  The  plate-like  morphology  of  clay  minerals  (e.g.. 
Fig.  2A)  contributes  to  the  persistence  of  these  terrigenous 
particles  in  water  columns  (Lai  and  Lerman  1975;  Davies- 
Colley  et  al.  2003).  Calcite  (not  a  noteworthy  component  in 
this  system),  in  contrast,  is  formed  autochthonously  during 
summertime  events  in  many  productive  hardwater  lacus¬ 
trine  systems  (Effler  1996;  Davies-Colley  et  al.  2003)  and 
has  been  linked  to  primary  production  (Hodell  et  al.  1998). 
This  mineral  can  become  dominant  with  respect  to  light 
scattering  during  these  events  (Peng  and  Effler  2005;  Peng 
et  al.  2007),  and  is  irregularly  important  in  three  of  the 
other  Laurention  Great  Lakes — Lakes  Ontario,  Erie,  and 
Michigan  (Strong  and  Eadie  1978). 

Uncertainties  in  SA X-Mie  estimates  of  bm  and  bhm — 
Sources  of  uncertainty  related  to  the  estimates  of  bm  and 
7>b  m  have  been  addressed  previously  (Peng  et  al.  2007).  A 
selective  review  of  that  material  is  included  here.  Because 
this  represents  the  first  effort  to  pursue  closure  for  bbp 
based  on  SAX-Mie  estimates  of  bbjns  particular  attention  is 
given  to  the  contrasting  situation  for  uncertainties  related 
to  bm  vs.  bh  m.  Issues  considered  here  are  limitations  of  the 
SAX  measurements  and  the  applicability  of  Mie  theory. 

The  morphological  characterizations  of  individual  par¬ 
ticles  by  SAX  are  robust.  The  precision  of  SAX  determi¬ 
nations  of  PAVm,  based  on  replicate  analyses,  is  about 
10%.  Only  modest  uncertainty  in  PA  (±5%),  and  related 
scattering  estimates,  is  attributable  to  the  protocol  adopted 
to  adjust  for  the  lying-flat  effect  (Peng  and  Effler  2007). 
There  is  only  minor  uncertainty  in  the  placement  of  SAX- 
characterized  minerogenic  particles  within  the  adopted 
chemical  classes.  However,  these  groups  are  not  definitive 
for  specific  (particularly  clay)  minerals,  thereby  introducing 
uncertainty  in  the  specification  of  the  m  value  (more 
specifically,  //).  The  most  different  of  the  common  clay 
minerals  in  this  regard  is  montmorillonite  that  has  an  n 
value  about  0.045  units  lower  than  most  of  the  others 
(Wozniak  and  Stramski  2004).  The  composition  of  the  ‘Si- 
rich’  and  the  unspecified  ‘Misc’  groups  (both  with  assumed 
m  values)  represents  another  minor  source  of  uncertainty, 
because  these  groups  only  made  modest  contributions  to 
PAVm  (Table  1).  Mie  theory  estimates  of  bm( 650)  have 
been  reported  (Peng  et  al.  2007)  to  be  insensitive  to  a  wide 
range  of  n  values  (broader  than  the  possible  deviation  for 
clays)  for  the  observed  PSD  patterns  (Fig.  3)  Increased 
sensitivity  is  predicted  for  Junge  PSD  cases.  In  contrast, 
estimates  of  6b>m  are  sensitive  to  modest  differences  in  n  for 


the  observed  PSD  pattern  (Peng  ct  al.  2007).  For  example, 
if  montmorillonite  was  dominant  instead  of  other  clay 
minerals  (as  assumed  here),  the  estimate  of  6bm  could  be 
biased  —20%  high  (see  fig.  9  in  Peng  et  al.  2007). 

Most  marine  (Volten  et  al.  1998;  Gordon  and  Du  2001) 
and  freshwater  (Peng  et  al.  2007;  this  study)  particles  are 
nonspherical,  and  thus  do  not  meet  the  sphericity  stipulation 
of  Mie  theory.  Deviations  from  sphericity  cause  shifts  in  the 
VSF  that  influence  light  scattered  in  the  backward  direction 
more  than  the  forward  (Mishchenko  and  Travis  1994). 
Depending  on  the  details  of  particle  shapes,  nonspherieity 
can  cause  either  enhancement  or  diminishment  of  the 
magnitude  of  bbp  relative  to  spherical  particles  (Bohren 
and  Huffman  1983;  Volten  et  al.  1998;  Gordon  and  Du 
2001).  The  effects  of  this  on  bp  are  genen  lly  less  (Holland 
and  Gagne  1970;  Jonasz  1987).  This,  together  with  the 
greater  sensitivity  of  Abp(A)  to  the  value  of  n  results  in  greater 
confidence  in  the  estimates  of  bm  that  those  of  /?b  m. 
However,  Liu  et  al.  (2003)  demonstrated  tl  at  for  nonspher¬ 
ical  quartz  particles,  the  deviation  of  the  M  e-modeled  phase 
function  from  that  of  measured  had  very  little  effect  on  the 
predicted  b^p.  More  complex  models  lor  nonspherical 
particles  (Mishchenko  and  Travis  1994)  vould  also  have 
substantial  uncertainty,  given  the  varia  bility  in  particle 
shapes  encountered  in  our  studies  (Table  1  Peng  and  Effler 
2007;  Peng  et  al.  2007).  Mie  theory  calci  lations  presently 
represent  the  most  viable  alternative  to  estimate  6ni  and  /;b#I11 
based  on  single  particle  optics  from  in  iividual  particle 
information  for  large  heterogeneous  part  cle  populations. 
The  applicability  of  Mie  theory  can  be  exp  xted  to  be  more 
tenuous  where  greater  deviations  fron  sphericity  are 
prevalent,  conditions  commonly  encounte  ed  for  biological 
particle  assemblages  such  as  pennate  diat(  ms  (Fig.  2B;  see 
ASP  values  for  diatoms  in  Table  1).  The  r  (bust  representa¬ 
tions  of  particle  shape  available  from  SAX  may  prove  to  be 
valuable  in  advancing  more  complex  light-  cattering  models 
for  minerogenic  particles. 

Pursuit  of  optical  closure— Two  other  jui  ces  of  uncer¬ 
tainty  need  to  be  considered  in  evaluati  ig  the  extent  of 
closure:  the  adequacy  of  the  empirical  bi  )-optical  models 
for  bQ  and  bbo,  and  the  representat  veness  of  bulk 
measurements  of  bp  and  bbp  and  collect  d  samples.  The 
models  (Eqs.  3  and  4)  developed  for  Case  oceanic  waters, 
driven  by  [Chi],  have  demonstrated  noteworthy  variability 
in  performance  for  those  targeted  ecosys  ems  (Morel  and 
Maritorena  2001).  A  major  contributor  to  this  variability  is 
the  imperfect  character  of  [Chi]  as  a  sur  ogate  of  optical 
properties  (Morel  and  Maritorena  2001).  The  dispersion  in 
the  />p-[Chl]  relationships  in  Case  1  \  aters  has  been 
described  as  “definitely  large”  (p.  7177;  Morel  and 
Maritorena  2001).  Substantially  greatci  uncertainty  in 
estimates  of  7>bi0  compared  to  that  of  bQ  s  acknowledged. 
The  exceedence  of  bbp  by  bb^0  estimates  at  ELS1  and  ELS2 
(Fig.  5C)  suggests  overestimation  of  the  organic  compo¬ 
nent  in  these  cases  by  the  bio-optical  nodel  (Eq.  4b). 
Divergent  /?bp-[Chl]  relationships  have  been  published,  but 
without  compelling  validation  (Morel  and  Maritorena 

2001).  Use  of  an  alternative  relationship  (Huot  et  al. 
2008)  developed  from  paired  Z?bp(650)  arJ  [Chi]  observa- 
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tions  in  Case  1  waters  caused  —40%  decreases  in 
predictions  of  bbo ,  but  resulted  in  similar  performance  in 
explaining  the  observed  variability  in  bulk  measurements 
with  the  two-component  model  (Fig.  6B). 

This  analysis  in  this  study  invokes  a  scattering  component 
for  organic  particles  in  Lake  Superior  that  corresponds  to 
Case  1  conditions.  The  range  of  [Chi]  in  the  lake  is  within 
that  covered  in  the  development  of  the  bio-optical  models. 
Other  covarying  scattering  particles  potentially  embedded 
within  these  relationships  include  organic  detritus  and 
minerogenic  particles.  The  assumptions  implicit  in  adoption 
of  Case  1  relationships  in  our  two-component  model  are  that 
any  embedded  detritus  is  representative  for  the  lake,  and  that 
any  embedded  minerogenic  component  is  insignificant 
compared  with  the  independent  SAX-Mie  estimates  of  bm 
and  bbm  (i.e.,  the  extent  of  ‘double  counting’  of  this 
component  is  negligible). 

The  variability  in  the  multiple  observations  of  bp  and  bbp 
within  the  depth  interval  of  measurements  was  relatively 
minor,  with  average  CVs  of  2.6%  and  5.5%,  respectively. 
The  uncertainty  in  the  bulk  measurement  of  bbpy  associated 
with  estimation  based  on  a  single  backward  angle 
measurement  (—10%;  Boss  and  Pegau  2001),  is  greater 
than  for  bp.  Whitmire  et  al.  (2007)  evaluated  multiple 
sources  of  uncertainty  for  bbp(/,)  and  presented  a  likely 
maximum  error  in  the  ratio  of  —20%,  a  magnitude  that  is 
modest  relative  to  the  observed  variability  (Fig.  5D). 
Finally,  the  widely  acknowledged  patchy  distribution  of 
particles  (Babin  et  al.  2003)  results  in  the  imperfect 
representation  of  bulk  conditions  and  by  water  samples 
collected  in  temporal  and  spatial  proximity. 

The  credibility  of  the  SAX-Mie  estimates  of  minerogenic 
scattering  in  Lake  Superior  is  supported  by  their  consis¬ 
tency  and  the  extent  of  closure  with  bulk  measurements 
through  the  two-component  approach  (Figs.  5B-D,  6). 
This  performance  for  such  a  dilute  system  and  the  addition 
of  closure  analyses  for  backscattering  (Figs.  5C,  6B) 
substantially  augments  earlier  efforts  to  test  the  credibility 
of  the  SAX-Mie  approach  in  other  lacustrine  systems 
(Peng  and  Effler  2007;  Peng  et  al.  2007).  The  weight  of 
continuing  successful  applications,  over  an  increasingly 
diverse  array  of  systems,  provides  the  strongest  case  for  the 
credibility  and  applicability  of  the  approach.  This  study  did 
not  prov  ide  a  rigorous  test  of  the  credibility  of  bm  estimates 
for  the  ELS  sites  because  of  the  minor  contributions  of  this 
component  to  bp.  However,  the  differences  in  this 
component  between  these  sites  and  those  in  the  embayment 
and  at  KW1  were  reasonably  well-supported  (Figs.  5B, 
6A).  The  greater  uncertainties  in  hbrn  estimates,  and 
contributing  components  of  related  closure  analyses  of 
bbo  estimates  and  bbp  measurements,  are  compensated  for 
by  the  stronger  signal  of  minerogenic  particles  for  /?bp 
compared  to  bp.  The  extent  of  closure  for  bbp  (Fig.  6B)  was 
actually  somewhat  better  than  that  for  bp  (Fig.  6A).  Much, 
if  not  all,  of  the  observed  shortcomings  in  closure  is  well 
within  the  accepted  variability  in  the  empirical  bio-optical 
models  for  Case  1  waters  adopted  here  for  the  organic 
components.  Moreover,  lhe  paradigms  describing  the 
dependence  of  bbp  on  the  relative  contributions  of  organic 
and  inorganic  particles  to  the  overall  particle  assemblage 


were  well-supported  by  the  results  from  tie  SAX  Mie 
approach  (Fig.  7C,D). 

The  demonstrated  extent  of  closure  has  implications 
with  respect  to  pathways  to  assess  scattering  conditions  in 
Lake  Superior  as  well  as  other  systems.  First,  reasonable 
first  approximations  of  the  organic  component  can  be 
made  from  widely  available  [Chi]  data.  Paired  observa¬ 
tions  of  PAVm  and  [Chi]  through  this  simple  two- 
component  model,  can  be  used  to  estimate  the  bulk 
scattering  metrics  of  bp  and  /?bp,  and  partition  these  into 
organic  and  multiple  minerogenic  components.  The  avail¬ 
ability  of  SAX  observations  is  presently  limited  However, 
these  findings  further  suggest  that  reasonable  fi rst  estimates 
of  bm  and  bb  rn  can  be  made  for  Case  2  waters  enriched  in 
minerogenic  particles  for  which  observations  of  bpy  6bp,  and 
[Chi]  are  available  (Loisel  et  al.  2007;  Whitmire  et  al.  2007), 
based  on  residual  calculations  {bbm  ~  bbp  -  bb%0). 
Conversely,  bQ  and  bbo  can  be  estimated  from  SAX  and 
bulk  bp  and  /?bp. 

Backscattering  ratio  at  650  nm  This  i  tio,  />bp,  is 
independent  of  particle  concentration,  with  \  hich  bp  and 
bbp  co-vary  linearly,  It  depends  instead  on  he  chemical 
composition,  PSD,  and  particle  shape(s)  of  the  assemblage. 
Thus,  bbp  represents  a  source  of  information  for  these 
attributes.  Measurements  of  bbp  have  been  use  !  to  describe 
the  variability  of  the  bulk  refractive  index  (;  )  of  oceanic 
particle  populations  (Twardowski  et  al.  2001  Boss  ct  al. 
2004;  Sullivan  et  al.  2005).  Low  n  values  a  rrespond  to 
dominance  by  phytoplankton,  whereas  high  '  allies  reflect 
dominance  by  minerogenic  particles  (Loisel  et  al.  2007; 
Whitmire  et  al.  2007).  In  addition  to  support!  ig  inferences 
of  particle  composition,  the  ratio  is  valuable  in  modeling 
the  underwater  light  field  (Mobley  et  al.  2002  . 

The  number  of  observations  of  bbp  for  this  j  tudy  is  small 
compared  with  surveys  reported  for  several  murine  systems 
(Loisel  et  al.  2007;  Whitmire  et  al.  2007;  Snyde  et  al.  2008). 
However,  our  observations  are  interesting  in  t  le  context  of 
the  strength  of  the  observed  relationships  with  [Chi] 
(Fig.  7A)  and  its  contribution  to  total  partic  ulate  matter 
([Chl]:rp;  Fig.  7B),  and  the  availability  »f  the  two- 
component  model  estimates  of  the  ratio  (1  ig.  7C)  and 
direct  minerogenic  scattering  (Fig.  7D).  T  le  two-fold 
difference  in  the  ratio  reported  here  was  sii  lilar  to  that 
observed  in  the  survey  of  several  Central  New  York 
systems  (Peng  et  al.  2007).  The  variability  witl  in  the  upper 
waters  of  Crater  Lake  was  apparently  less,  w  th  only  17% 
difference  between  the  10th  and  the  90th  percentiles  of  the 
observations  (Whitmire  et  al.  2007).  Much  br  nider  ranges 
(e.g.,  >  10-fold  differences)  in  bbp  were  repor  ed  in  larger 
marine  surveys  (Boss  et  al.  2004;  Loisel  ;t  al.  2007; 
Whitmire  et  al.  2007).  However,  the  central  tendency  of 
the  Lake  Superior  observations  was  similar  to  the  other 
studies;  for  example,  the  mean  was  0.0108,  cc  npared  with 
0.010  reported  by  Whitmire  et  al.  (2007)  and  0  0138 
(geometric  mean)  by  Loisel  et  al.  (2007). 

The  performance  of  the  power-law  fit  R2  0.65, 
Fig.  7A)  describing  the  dependence  of  /),>p  on  [Chi] 

compares  favorably  with  those  reported  elset  here  (Twar¬ 

dowski  et  al.  2001;  Sullivan  et  al.  2005;  Wl  itmirc  et  al. 
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2007).  However,  this  model  is  mechanistically  weak,  and  in 
this  case  inappropriate,  because  it  fails  to  accommodate  the 
documented  spatial  differences  in  minerogenic  scattering 
(Fig.  5C)  that  did  not  covary  with  [Chi]  (Fig.  5A).  The 
dependence  of  bbp(X)  on  the  ratio  [Chi] :  cp(X)  is  more 
mechanistically  sound,  because  this  ratio  serves  as  a  proxy 
for  the  relative  proportion  of  phytoplankton  within  the 
overall  particle  assemblage  (Boss  et  al.  2004;  Loisel  et  al. 
2007),  and  it  was  stronger  ( R 2  =  0.88;  Fig.  7B).  Decreases 
of  bbp  reflect  increases  in  the  contribution  of  phytoplankton 
to  the  light-scattering  particle  assemblage  (Boss  et  al.  2004; 
Loisel  et  al.  2007).  The  slope  of  the  Lake  Superior 
relationship  (-0.0044)  is  similar  to  that  (-0.0066)  reported 
by  Boss  et  al.  (2004)  between  observed  6bp(532)  and 
[Chi] :  cp(660).  The  relationship  remained  strong  for  the 
backscattenng  ratio  based  on  the  two-component  model 
(Fig.  1C). 

The  greater  strength  of  the  observed  bbp  vs.  [Chi] ;  cp 
relationship  ( R 2  =  0.88)  for  the  single  cruise  compared  to 
larger  marine  surveys  is  noteworthy.  For  example,  the 
relationships  reported  by  Boss  et  al.  (2004)  and  Loisel  et  al. 
(2007)  explained  56%  and  26%  of  the  observed  bbp 
variability.  This  suggests  greater  uniformity  in  the 
monitored  light-scattering  features  of  particles  in  Lake 
Superior  for  this  single  cruise.  A  contributing  factor  is  the 
similar  light-scattering  characteristics  of  the  minerogenic 
particle  populations  other  than  A  (Figs.  3,  4,  5E;  Table  1). 
The  strength  of  the  observed  bbp  -  [Chl]:cp  relationship 
further  suggests  a  relatively  uniform  cellular  content  of 
chlorophyll  a  for  the  narrow  [Chi]  range  observed  here. 
Greater  variability  in  this  relationship  (Fig.  7B)  may  arise 
from  year-round  monitoring  as  the  cellular  content  of  this 
pigment  (Reynolds  2006)  shifts,  even  if  the  light-scattering 
attributes  of  the  minerogenic  particles  remain  unchanged. 
However,  it  remains  that  over  this  brief  interval  the  ratio 
[Chi] :  cp  was  a  strong  predictor  of  observed  bbp.  The  Lake 
Superior  results  suggest  that  the  relationship  may  be 
substantially  stronger  within  individual  systems  than  those 
observed  for  broad  cross-sectional  studies  (Loisel  et  al. 
2007;  Whitmire  et  al.  2007). 

A  corollary  interpretation  of  the  observed  relationship 
between  />bp  and  [Chi] :  cp  (Fig  7B)  described  elsewhere  has 
been  that  the  contribution  of  minerogenic  particles  to  the 
light-scattering  particle  assemblage  increases  as  bbp  increas¬ 
es  (Twardowski  et  al.  2001;  Boss  et  al.  2004;  Loisel  et  al. 
2007).  Although  this  position  has  been  supported  by 
corresponding  estimates  of  bulk  refractive  index,  it  has 
not  been  rigorously  supported  by  paired  direct  measure¬ 
ments  of  the  inorganic  particle  population.  The  estimates  of 
bm  provide  the  opportunity  to  directly  test  that  interpre¬ 
tation  for  this  study  Indeed,  a  strong  (R2  —  0.67)  positive 
relationship  prevailed  between  bbp  and  the  ratio  bm  :cp  for 
this  study  (Fig.  6D),  directly  supporting  the  position  that 
bbp  increases  as  the  fraction  of  light  scattering  associated 
with  minerogenic  particles  increases.  A  similarly  strong 
relationship  was  observed  when  the  bm  estimates  in  the 
ratio  were  replaced  by  the  PAVm  measurements.  The 
calculated  bbp  values  for  hypothetical  cases  of  strictly 
minerogenic  or  phytoplankton  populations  (Fig.  5E)  can 
be  considered  as  bounding  conditions  for  this  system.  The 


observed  intermediate  values  (Fig.  5D)  represented  the 
outcomes  of  mixtures  with  varying  contributions  from  the 
minerogenic  and  organic  (phytoplankton)  particle  popula¬ 
tions. 

Implications -  The  integrated  SAX-Mie  approach  offers 
great  promise  for  enhancing  the  understanding  of  the  role 
of  minerogenic  particles  in  influencing  the  underwater  light 
field  and  remote-sensing  signals.  Applications  of  this 
approach  advance  implementation  of  the  reductionist 
strategy  for  better  understanding  of  the  real  complex 
optical  regimes  through  quantitative  partitioning  of  the 
effects  of  organic  vs.  multiple  minerogenic  constituents. 
Inclusion  of  SAX  characterizations  m  the  suite  of 
measurements  of  optics  monitoring  programs  would 
provide  independent  measurements  of  PSD  and  composi¬ 
tion  for  multiple  constituents  within  minerogenic  particle 
populations,  thereby  avoiding  the  need  for  arbitrary 
specifications  in  related  modeling  efforts  This  represents 
an  invaluable  reduction  in  the  'degrees  o  freedom’  of  the 
modeling  process  by  providing  indepen  lent  key  inputs, 
widely  sought  by  the  modeling  community  to  improve  the 
credibility  of  such  frameworks  (Chakra  1997).  SAX 
information  has  potential  value  for  a  wid  range  of  issues, 
such  as  to  address  origins  of  backscatt ‘ring  in  oceanic 
waters  (Stramski  et  al.  2004a),  and  to  ai  i  in  establishing 
management  targets  and  reasonable  <  xpectations  lor 
programs  intended  to  increase  clarity  in  inland  waters 
(Effler  et  al.  2008). 

This  study  of  Lake  Superior  has  enhai  ced  the  case  for 
the  credibility  of  the  SAX-Mie  approach  o  independently 
quantify  bm,  and  has  provided  a  valuable  md  encouraging 
first  test  of  the  capability  for  6b,m-  Lake  Superior  has 
represented  a  challenging  test  for  inland  water  application 
because  of  the  relatively  dilute  nature  of  this  system, 
thereby  providing  optimism  for  successfu  application  for 
more  spatially  and  temporally  robust  descriptions  of  this 
lake  and  for  other  dilute  systems,  including  the  other 
Laurentian  Great  Lakes  and  oceanic  w  ters.  Moreover, 
increased  support  for  bbp  as  a  surroga  e  metric  of  the 
contribution  of  minerogenic  particles  to  <  verall  scattering 
has  been  provided. 

Much  broader  applications  of  the  SA)  Mie  approach, 
in  concert  with  the  bulk  measurements  of  bp  and  bbp 
included  here,  are  recommended  to  (1)  further  test  the 
consistency  and  related  extent  of  closure  lor  these  forward 
estimates  of  light  scattering,  (2)  demonsti  ite  the  utility  of 
the  increased  partitioning  capability  of  the  minerogenic 
component,  and  (3)  advance  the  testing  of  Jbp  as  a  metric  of 
the  relative  contributions  of  mineroge  lie  particles  to 
overall  scattering.  Increased  testing  is  parti  ;ularly  necessary 
for  SAX-Mie-based  estimates  of  bbm  to  alliance  the  case 
for  the  credibility  of  these  estimates  a  id  to  represent 
associated  effects  on  the  remote-sensing  s  gnal. 
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